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Abstract 

The rare decays Bs — ?• /i"^/^^ and B — t- K^*^ fj,~ are sensitive to new particles and 
couplings via their interferences with the standard model contributions. Recently, the 
upper bound on B{Bs — )• has been improved significantly by the CMS, LHCb, CDF, 

and D0 experiments. Combining with the measurements of B{B — t- K^*^ fj,~^ fi~), we derive 
constraints on the relevant parameters of minimal supersymmetic standard model with 
and without R-parity, and examine their contributions to the dimuon forward-backward 
asymmetry in — t- K*fi^fj,~ decay. We find that (i) the contribution of R-parity violating 
coupling products X2i2^2i3 "^^^ t° squark exchange is comparable with the theoretical 
uncertainties in i? — )• K^^ decay, but still could be significant in B — )• K*^^fj,~ decay 
and could account for the forward-backward asymmetry in all dimuon invariant mass 
regions; (ii) the constrained mass insertion ((^j^^)23 could have significant contribution to 
dApsiB —7- K* ^'^ fj,'~)/ds, and such effects are favored by thr recent results of the Belle, 
CDF, and LHCb experiments. 
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1 Introduction 



Recently, using the 7/6 ^ data set, the CDF Collaboration at the Fermilab Tevatron has ob- 
served an excess of Bg candidates which is compatible with 



(1) 



and provided the corresponding upper limit of B{Bs — )■ fi^fi ) < 4.0 x 10 ^ at 95% confidence 
level (CL). 

At the same time, searches for Bg — )■ fi~^fi~ have also been made by the CMS and LHCb 



Collaborations 



^, respectively, at the Large Hadron Collider at CERN. The combined results 



of the searches by the CMS and LHCb Collaborations in the upper limits jsj] are 

B{Bs < 1-08 X 10^^ at 95% CL, (2) 

B{Bs < 0.90 X 10"^ at 90% CL, (3) 

which have improved the previous upper bounds ^\ significantly. 

Bs — )■ P-~^fJ'~ decay is a known sensitive probe to the presence of new physics (NP). In the 
standard model (SM), it occurs via penguin or box diagrams and is strongly helicity suppressed. 
Its SM prediction is (3.2 ± 0.2) x 10"^ Generally, NP could enhance the B^ -> 
decay rate very much, and thus the upper bound of B{Bs — )■ fi~^fi~) is taken as a strong 
constraint when a NP model is discussed. As a cross-check, one usually needs to investigate the 
semileptonic rare decays B — Kfi~^fi~ and B — j- K*fi^fi~ which are also governed by the flavor 
changing neutral current transition b — sfi'^fi" but not helicity suppressed. Many observables 
of i? — )■ K^*^fi'^n~ have been observed by several experiments: BABAR 3, Belle 3, CDF 



10|, and LHCb Uj. As many of them agree with the SM predictions within their error bars. 



however, the dimuon forward-backward asymmetry of -B — ?■ K*fi'^fi at the low region o: 
dimuon invariant mass is not consistently measured by Belle js], CDF 10 1, and LHCb 11|. 



the 



Any NP that alters B{Bs — ?■ would necessarily alter observables in i? — > K^*^ii^ii~ 

decays; examples of the latter are the differential branching ratio and forward-backward asym- 
metry. The NP effects in the b — )■ sfi~^fi~ flavor changing neutral current transition have been 



extensively inves 



tigated, for instance, in Refs. 12N21|. In this paper, following closely the 



analysis of Ref. 22j, we will update the constraints on the R-parity violating (RPV) minimal 



23L 24 1 ■ Using a com- 
5| as well as the 



supersymmetric standard model (MSSM) in light of the new experimental data on Bs — )■ /i^/i 
and B K^*^ jX^H . Additionally, we will extend our analysis to the R-parity conserving 
(RPC) MSSM scenario with the mass insertion (MI) approximation 
bination of the limits of from CDF, LHCb and CMS 

experimental bounds of B{B — )■ K'^*^ fi'^ 25|, we will obtain the new limits on the relevant 
supersymmetric coupling parameters. Then we will use the constrained parameter spaces to 
examine their effects on some observables in these decays, especially dApsiB — )■ K*fi~^fi~)/ds. 

The paper is arranged as follows. In Sec. 2, we present a very brief theoretical introduction 
to Bs — and B — t- K'^*^ jX^ fj.'^ processes. In Sec. 3, we deal with the numerical results. 
We display the constraints implied by the new experimental data on the RPV and RPC MSSM 
parameter spaces and discuss the implications for the B — K^*^ fX^ fj," invariant mass spectra 
and forward-backward asymmetries. Section 4 contains our conclusion. 



2 The theoretical framework for Bg — > /i^/i and B 
decays 



2.1 The leptonic decay Bg 



The branching ratio for Bg — fi^fi can be written as 19|, |26| 
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Fb 

The SM result for the branching ratio may be obtained from Eq. 
Cs = Cp = Cp = and 

Gpae 



C, 



VtbV*Y{x,). 



(4) 



(5) 



by setting Ca = Cs 



(6) 
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In the MSSM without R-parity, the branching ratio may be obtained by setting 22 

C'a = 



C' 



2 ' 



c' 



_ Ai22A-2 3 
^i22"^t32 



(7) 



In the MSSM with R-parity, the branching ratio can obtained by using the expressions 
Cs, Cs, Cp and Cp can be found in Ref. [19|; and Ca = in this case. 



2.2 The semileptonic decays B — > K^*^^!^ 

In the SM, the double differential decay branching ratios and ^^^7^ for the decays B — )■ 



Kn'^H and B — )• K*fi~^fi , respectively, may be written as 



27| 
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where p = Pb + Pk(*)^ s = q'^, and q = p+ + p^ ( p+ the four-momenta of the muons), and the 



auxiliary functions A — H can be found in Ref. 



271 ]. The hat denotes normalization in terms 



of the B- meson mass, mp, e.g., s = s/m^, rhq = mq/mp. 
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In the MSSM without R-parity, the double differential decay branching ratios including the 
squark exchange contributions could be gotten from Eqs. I^IM) by the replacements |22| 
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where W = -^V*Vt!,7nB. 

The sneutrino exchange contributions are summarized as 



(10) 
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In the MSSM with R-parity, all the effects arise from the RPC Mis contributing to C7, C^^^ , Cio, 
and they are 
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where n = 1 for decay B Kfi'^fi as well as for the terms related to the form factors V 
and Ti in i? — i- K*fi'^fi'' decay, n = —1 for the terms related to the form factors Aq, Ai, A2, T2 
and T3 in 5 ^ JC/i+zi" deca^C:^'"^'^^ (Cg^^^)^*'^^''^^^ C'fo"^'*^^ Ci/^', (C^^^*'' and C[^' 



have been estimated in Refs. 



28 



obtained from Eqs. ([319]) by the following replacements [l^, H]- 



30|. The results for and S^* including MI effects can be 
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From the total double differential branching ratios, we can get the dimuon forward-backward 
asymmetries 



+1 d'^BiB^RMfi+fj.") 



dsdcosd 
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(16) 



3 Numerical results and analyses 



We will present our numerical results and analysis in this section. When we study the effects 
due to MSSM with and without R-parity, we consider only one new coupling at one time, 
neglecting the interferences between different new couplings, but keeping their interferences 
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with the SM amphtude. The input parameters are collected in the Appendix, and the fol 
experimental data will be used to constrain parameters of the relevant new couplings 

B{Bs < 1-08 X 10^^ (at 95% CL), 

B{B Kfi-^fi") = (0.48 ± 0.06) x 10~^ 
B{B = (1-15 ± 0.15) X 10"^ 



owmg 

urn: 



(17) 



To be conservative, we use the input parameters varied randomly within la variance and the 
experimental bounds at 95% CL. We do not impose the experimental bounds from dApBiB — )■ 
K* fi~)/ds and leave it as predictions of the restricted parameter spaces of the two NP 



scenarios, and compare them with the experimental results in Refs. j9-lll|. 



3.1 RPV MSSM effects 

First, we will consider the RPV effects and further constrain the relevant RPV couplings from 
the new experimental data of B{Bs — )■ /i^yU~) and B{B — > K^*^ ji'^ given in Eq. f[T7|) . As 
given in Sec. El there are three RPV coupling products, which are )^2i2'^2a due to squark 
exchange as well as \i22^'i2z and A*22A'j32 due to sneutrino exchange, relevant to Bg 
and B — > Ji decays. 

Our new bounds for three RPV coupling products from the 95% CL experimental data 
are demonstrated in Fig. [H And the upper limits for the relevant RPV coupling products 
by B{B — )■ K^*^ and B{Bs — t- /i'^'yu") are summarized in Tabled] For comparison, our 
previous bounds on these quadric coupling products are also listed. From Fig. [T] and Table 
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Figure 1: The allowed RPV parameter spaces with 500 GeV sfermions, and the RPV weak 
phase {(pRpv) is given in degree. 
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Table 1: Bounds for the relevant RPV coupling products hj B K^*^ fj,'^ fj, and Bg — )• fj,^ fj, 
decays for 500 GeV sfermions, and previous bounds are listed for comparison. 



Couplings 


Bounds 


Previous bounds [22] 


1 \' V* 1 

Aj22A'j32 

W* V 1 

I^i22^i23l 


< 8.2 X 10-^ 

< 2.0 X 10^^ 

< 2.0 X 10^^ 


< 11.5 X 10"^ 

< 4.5 X 10"^ 

< 4.3 X 10-^ 



HI one can find that all three RPV coupling products are restricted, and the upper limits of 
|Ai22A'j32| and |A*22A^23l are improved by about a factor of 2 by the new experimental data. 
Notice that we assume the masses of sfermions are 500 GeV. For other values of the sfermion 
masses, the bounds on the couplings in this paper can be easily obtained by scaling them by 
factor of p = (soSfcy)'- 

Now we will analyze the constrained RPV effects on B{Bs — )■ yU^/i"). The sensitivities of 
B{Bs — )■ fJ.'^fJ.') to the constrained RPV couplings are shown in Fig. [2J The limits of the 
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Figure 2: The constrained RPV coupling effects on B{Bs — )■ /i^/i ). The olive (violet) 
horizontal dotted (solid) lines denote the limits of the 95% CL measurements (SM predictions). 



8 



measurements at 95% CL and the SM predictions with la theoretical uncertainties are also 
displayed in Fig. |2] for convenient comparison. Figs. [2] (a) and (b) show the constrained 
effects of the modulus and weak phase of t-channel squark exchange coupling A2j2-^2i35 respec- 
tively. As shown in Figs. [2] (a-b), with the contribution of Ag^gA'g^g included, B{Bs /U"*"/^") is 
lower than its experimental upper limit 5|. Besides the constraints from B{Bs — )■ K^*^ 
^2i2^2i3 coupling is uot further constrained by the new experimental upper limit from CMS 
and LHCb since its contribution to B{Bs — is suppressed by w?^/m\. Additionally, 
the allowed parameter space of A',,A& would be excluded tt the 68% CL experiuteutal detet- 
mination B{Bs fi+fx') = il-8tl:l) x 10"^ y by the CDF Collaboration were taken as a 
constraint. Two s-channel sneutrino exchange contributions to B{Bs — )■ fJ'^fi") are very similar 
to each other. We would take the Xi22K*23 contribution as an example, which is shown by Figs. 
|2](c-d). We can see that B{Bs —J- is sensitive to both the modulus and phase of Ai22A'j23, 

and B{Bs — ?■ not only could be increased but also could be decreased by the presence of 

At22A^23 coupling. Generally, the Ai22A^23 coupling could alter B{Bs — )■ fi^fi~) significantly since 
its contribution is not helicity sirppressed by m^/m^. Thus, the constraint on Aj22A^23 is due 
to the bound of B{Bs — > [5|. 

Then we turn to analyzing the constrained RPV effects in B ^ K^*^ /i^/i" decays. Using the 
new constrained parameter spaces shown in Fig. [T], we will give the RPV effects on the dimuon 
invariant mass spectra and the forward-backward asymmetries of i? — Ji decays. 

In Fig. [31 we present correlations between the dimuon invariant mass spectra as well as 
the dimuon forward-backward asymmetries and the parameter spaces of A2j3A2*2 by the two- 
dimensional scatter plots. The dimuon invariant mass distribution and the dimuon forward- 
backward asymmetry are given with vector meson dominance contribution excluded in terms of 
dB/ds and dApB/ds, and included in dB' /ds and dA'p^/ds, respectively. In Fig. [3l the magenta 
"x" denotes the SM prediction within la error ranges of the input parameters, olive solid line 
denotes the central value of the SM prediction, and blue "|" denotes the RPV supersymmetry 
(SUSY) prediction including A2j2A2*t3 coupling within la error ranges of the input parameters. 
The theoretical uncertainties of the SM predictions of dB{B — )■ K^*^ fi~)/ds are quite large; 
nevertheless the theoretical uncertainties are canceled a lot in dApsiB — ?■ K* fi^ fi'')/ds. 

The RPV effects on dA'p^lB — t- K*fi~^fi~)/ds are shown in Fig. [3] (f). This observable 
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Figure 3: The effects of RPV coupling A2j2A2*3 due to the squark exchange in 5 — > K^*^n'^fi~ 
decays. The w denotes s, magenta "x" denotes the SM prediction within la error ranges of 
the input parameters, ohve sohd hne denotes the central value of the SM prediction, and blue 
"I" denotes the SUSY prediction. The same goes for Figs. lU [HI andO 

las been measured as a function of the dimuon invariant mass square by BABAR 8|, Belle 
9], CDF 10], and LHCb [ll|, and the current situation is specially exemplified in Fig. |H As 
shown in Fig. HI the fitted dA!p^{B — K*fj,~^fj,~)/ds from Belle is generally higher than the SM 
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Figure 4: Afb{B K*ii^n~) including RPV coupling X'2i2^2i?. versus the 95% CL data: CDF 
(blue dotted line), Belle (purple solid line), and LHCb (red dashed line). 

expectation in whole bins, and the CDF fitted result is consistent with the SM prediction 
in some bins and it is higher than the SM prediction in some other bins; nevertheless 
the LHCb fitted result, which is the most precise to data, is in good agreement with the SM 
prediction. Especially, in the region of < s < 0.072 (i.e., GeV^ < < 2 GeV^), the 
Belle measurement favors a positive value which is not confirmed by CDF and LHCb, whereas 
the sign of the SM prediction for dA!pg{B — )■ K*fi~^fi~)/ds is negative. One could find that 
the constrained RPV coupling A2j3A2*2 still could accommodate dApBiB — t- K*fi~^fi^)/ds from 
Belle, CDF, and LHCb at all s regions. 

As for the s-channel sneutrino exchange couplings Xi22K*23 -^122 -^'132; constraints 
from B{B — )■ fi^fi') are rather restrictive. The Aj22A^23 coupling effects in i? — )■ K^*^ii'^^~ 
are displayed in Fig. [Sj we see that Aj22A^23 coupling has negligible contribution to dB{B — )• 
fi+fi )/ds, and the differences between the SUSY prediction and the SM ones are due to 
the 95% CL experimental constraints. Nevertheless, constrained Aj22A'j23 coupling has some 
effects on dApsiB — )• K* n~) /ds. \i22K32 coupling effects in i? — > K^*^n'^n~ are similar to 
Aj22A^23 effects; thus we will not show them again. 



11 




Figure 5: The effects of RPV coupling \i22^'i2-A to the sneutrino exchange in B 
3.2 RPC MI effects 

Now we study RPC MI effects in B, n+jj' and B I&^n+jj- decays in the MSSM with 
large tan/3. The eight kinds of Mis {6^3)23 with {A,B) = {L,R) contribute to B ^ K^*^fi^fi~ 
decays, but only three kinds of Mis {S^l)23, {^11)23, and (5^^2)23 contribute to Bs —J- fi'^fi' decay. 
We will only consider the contributions of (5^^)23, (^ll)23, and (5^/j)23 Mis to Bg — yU^/i" and 
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Figure 7: The constrained MI effects on B{Bs — )■ fi^fi ). The ohve (violet) horizontal dotted 
(solid) lines denote the limits of the 95% CL measurements (SM predictions with la error bar). 

predictions at |(5^'^)23| = and the SM predictions come from contributions in the MSSM with 
the CKM matrix as the only source of flavor violation. 

Then we analyze the constrained (5^^)23, {^11)23^ and (5^r)23 MI effects in 5 — )■ K^*^fi^fi~ 
decays. Using the constrained parameter spaces shown in Fig. [6l we will give the MSSM predic- 
tions to the dimuon invariant mass spectra of the decay width and the dimuon forward-backward 
asymmetries of 5 — K^*^jj~^jj~ decays in the MI approximation. Besides the MI contributions, 
the SUSY predictions also include the contributions that come from graphs including SUSY 
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Figure 8: The constrained non-MI effects in 5 — ?■ K^*"* jX^ jj, decays. 



Higgs bosons and sparticles in the hmit in which we neglect all the MI contributions, which are 
called non-MI contributions, and the non-MI SUSY effects are shown in Fig. |8l From Figs. [8] 
(a-b), we can see that dB{B Kfi~^fi~)/ds could be slightly suppressed at all s regions by the 
non-MI SUSY couphngs. As shown in Figs. IH](c-d), dB{B K*fi'^fi~)/ds could be decreased 
a lot at the middle s region by these couplings. Figs. M (e-f) show us that the non-MI SUSY 
couplings could slightly suppress dApB^B — )■ K* fi~^fi~)/ds at the middle s region. 
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The constrained (5^^)23 and (5^/2)23 Mis have no obvious effects in i? — )■ K{*)fi~^fi decays. 
(5El)23 mi contributions to 5 ^ K^*^ fi'^fi are presented in Fig. |9l Note that the SUSY 
predictions in Fig. 8 also include the non-MI contributions shown in Fig. |8l As shown in Figs. 
I2](a-b), the constrained {S'^i)23 MI has no obvious effects on dB{B Kfi^ fi^)/ds, which could 
be slightly suppressed at all s regions by only non-MI effects. On the other hand, its contribution 
to i? — K*fi~^fi~ could be significant, as shown in Figs. [9] (c-f), when theoretical uncertainties 
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Figure 9: The constrained {5'^l)23 MI effects in 5 — )■ K^*^fi'^jj, decays. 
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are considered. It is of interest to note that the contribution to dApB^B — )■ K*jj,^_u )/ds is 
favored by the current experimental measurements from Belle, CDF, and LHCb 9l-lll|. 



4 Conclusions 

Motivated by the recent searches of B{Bs — )■ yu'^/i") by the CDF, LHCb, and CMS Collab- 
orations, we have studied Bs — t- and B K^*^jj,'^fi~ decays in the MSSM with and 
without R-parity. In the MSSM without R-parity, we have found that the bounds of sneutrino 
exchange RPV couplings are significantly improved by the present new measurements. The 
further constrained RPV coupling due to t-channel squark exchange still has significant effects 
in i? — K^*^fj.~^fx~ decays, and the current measurements of dApsiB — )■ K* / ds could 
be accommodated by the squark exchange coupling. The further constrained couplings due to 
s-channel sneutrino exchange could have large effects in B^ — )■ fi~^fi~, but have negligible effects 
in B ^ K^*^fi'^fi~ decays. 

In the MSSM with R-parity, three MI parameters {S^l)23, (^ll)23, and (5^/j)23 suffer the 
combined constraints from the present data of B{Bs — fJi^ fi") and B{B — > K^*^ The 
constrained {6^1)23 MI could give large contributions to dApB^B — K* ji^ ji") / ds at all s 
regions in favor of the current experimental measurements from Belle, CDF, and LHCb. The 
constrained {Sfi j^i) 23 Mis have ignorable effects on the observables of i? — A'^*)/i"''yU~ decays. 
dApB^B — )■ K*ii^ji~)/ds could be slightly decreased at the middle s region by the SUSY 
contributions which come from graphs including SUSY Higgs bosons and sparticles in the limit 
in which we neglect all the MI contributions. 

In the immediate future, the LHC is expected to become sensitive to B{Bs jX^n'). 
Accurate measurements of the Bs — ?■ and B — )■ K^*^fi~^fi~ decays could further shrink or 

reveal the parameter spaces of MSSM with and without R-parity. 
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Appendix: Input parameters 



The input parameters are summarized in Table [2j For the RFC MI effects, we tak;e the five 
free parameters rriQ = 450 GeV,mi/2 = 780 GeV,AQ = —1110, sign(/i) > and tan/3 = 41 



from Ref. 



3l| . All other MSSM parameters are then determined according to the constrained 



MSSM scenario as implemented in the program package SUSFECT 3J]. For the form factors 



mvo 



35 



ving the B — )■ K^*^ transitions, we will use the recent light-cone QCD sum rules results 
which are renewed with radiative corrections to the leading twist wave functions and 



SU(3) breaking effects. For the dependence of the form factors, they can be parameterized 
in terms of simple formulas with two or three parameters. The expression can be found in Refs. 

In our numerical data analysis, the uncertainties induced by F{0) are also considered. 



35 



Table 2: Default values of the input parameters. 



TUB, = 5.370 GeV, m^, = 5.279 GeV, niB^ = 5.279 GeV, mw = 80.425 GeV, 

mK± = 0.494 GeV, tuko = 0.498 GeV, mK*± = 0.892 GeV, rriK^o = 0.896 GeV, 

m^K) = (4.191°:^^) GeV, m,(2GeV^) = (0.1001°;°^°) GeV, 

m„(2Ge1/) = 0.0017 ~ 0.0031 GeV, md{2GeV) = 0.0041 ~ 0.0057 GeV, 

me = 0.511 X 10-3 GeV, = 0.106 GeV, mt,poie = 172.9 ± 1.1 GeV [25] 



r^^ = (1.466 ±0.059) ps, r^^ = (1.530 ± 0.009) ps, r^^ = (1.638 ± 0.011) ps. 


[25] 


\Vtb\ ^ 0.99910, \Vts\ = 0.0416li[];S. 


[25] 


sin^^VK = 0.22306, ae = 1/137. 




fn, = 0.230 ± 0.030 GeV 


m 
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